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Using exact diagonalization in the spherical geometry, we investigate systems of localized quasi- 
holes at f = 5/2 for interactions interpolating between the pure Coulomb and the three-body 
interaction for which the Moore- Read state is the exact ground state. We show that the charge e/4 
quasihole can be easily locahzed by means of a (5-function pinning potential. Using a tuned smooth 
pinning potential, the quasihole radius can be limited to approximately three magnetic length units. 
For systems of two quasiholes, adiabatic continuity between the Moore- Read and the Coulomb limit 
holds for the ground state, while for four quasiholes, the lowest two energy states exhibit adiabatic 
continuity. This implies the existence of a Majorana fermion for pure Coulomb interaction. We also 
present preliminary results in the Coulomb limit for braiding in systems containing four quasiholes, 
with up to 14 electrons, diagonalizing in the full spin-polarized sector of the second Landau-level 
Hilbert space. 
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I. INTRODUCTION 

The possibility of realizing non-Abelian braiding 
statisticgH in a condensed-matter system has generated 
a great deal of interest both for theoretical and experi- 
mental studies. One such system that may be closest to 
experimental realization is the fractional quantum Hall 
(FQH) state at filling fraction ly — 5/2, first observed by 
Willett et aZ.El Whether the experimentally observed 5 /2 
state indeed has charged excitations with non-Abelian 
braiding statistics is still unknown, altough the discovery 
of a neutral current in experiment may make us hopeful.'^ 
For a review, we refer the reader to Nayak et al.^ and to 
Stern.lS 

In an earlier work® we presented numerical evidence 
that the FQH ground state (GS) at = 5/2 for Coulomb 
interaction, when spin polarized, is in the same univer- 
sality class as the non-Abelian Moore-Read (MR) Pfaf- 
fian state. ^ Calculating the energy spectrum of the few 
lowest-energy states by exact diagonalization for electron 
interactions interpolating between the Coulomb interac- 
tion Vc and the three-body interaction V^b, for which the 
MR state is the unique GS, we found for all examined 
system sizes adiabatic continuity (AC) of the GS and no 
sign of a decrease of the gap between the two limits Vc 
and Vsf). We concluded that AC can be expected even in 
the thermodynamic limit. In addition, we drew a phase 
diagram in the two-body interaction space (in the vicin- 
ity of the Coulomb interaction), showing that the FQH 
gapped phase coincides with the MR phase. 

Previous theoretical worlPltiSl gave indications for AC 
for the spin-polarized ly — 5/2 GS between the MR and 
the Coulomb limit. In the disk geometry Xin Wan et al.^ 
were also able to localize a single quasihole (QH) at the 
center of the disk showing that, in a certain range of QH 
pinning and system confinement potential strength, the 
lowest-lying state belongs to the total angular momentum 
value that corresponds to a MR e/4 charged QH, for 



electron interactions interpolating between V^b and Ve- 
in spite of the AC of the GS, the question whether 
the elementary charged excitations preserve their non- 
Abelian properties in going from the MR to the Coulomb 
limit still needs to be examined: This is a very important 
issue, because braiding of non-Abelian quasiparticles has 
been proposed for topological quantum computation.'^'ii' 
At 1/ = 5/2, one expects that QHs can have either charge 
e/4 with the possibility of having non-Abelian braiding 
statistics or e/2 with Abelian fractional statistics. In- 
deed, the results by Toke et al\^ have cast doubt on the 
existence of localized QHs with charge e/4. 

Here, we study fully spin-polarized systems with local- 
ized QHs in the spherical geometry. This polarization 
choice is motivated by theoretical investigations concern- 
ing this issue: One of us showed that the GS of the 
disorder-free FQH state at v = 5/2 is spin polarized 
even for vanishing Zeeman energyp^ recent theoretical 
worlfSmi confirmed this result. However, from the exper- 
imental point of view, the situation is less cle ar. T rans- 
port experiments with variable electron densitj^^^^ or in 
a tilted magnetic field at high density-* have been inter- 
preted in controversial ways: Varying the electron den- 
sity by more than a factor 2,1^ the excitation gap varies 
smoothly with no break in slope or discontinuity, indicat- 
ing that neither GS nor excitations change their character 
while the Zeeman energy varies by an amount that ex- 
ceeds the measured energy gap by a factor of 5 — 10. That 
an unpolarized state would survive under this condition 
appears unlikely. Also, in high-density samples, the dis- 
appearence of the gapped phase in a tilted magnetic field 
cannot be explained by a Zeeman energy term with a 
GaAs bulk g factor .1^^ Nevertheless, arguments have been 
presented that these experimental results might be con- 
sistent with a spin-unpolarized 5/2 state.'^^' Also first 
results from direct spin p olarization measurements, us- 
ing optical methods ,'221211 appear to be consistent with a 
spin-unpolarized state. As the experimental samples are 
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likely to be inhomogeneous and consist of compressible 
regions separated by percolating incompressible filaments 
at ly ^ 5/2,^^ it is not clear that the optical data relate 
solely to the filaments or to the compressible regions!^ 
The separate determination of the polarization state of 
the incompressible parts has not yet been possible. In 
addition, charged impurities may also induce skyrmions 
in the incompressible domains, which will decrease 
the total spin polarization. We also note that, once ther- 
mally assisted tunneling is taken into account ,122] the gap 
obtained from analyzing the dissipative conductance can 
be as consistent with expectations for a spin-polarized 
system at v = 5/2 (Refs. 17, 24 and 25) as for the hier- 
archy state at v — A/9. 

In our investigations at first we show that charge e/4 
QHs can be localized by a (5-function potential and, vary- 
ing the form of the localization potential away from the 
simple (5-function, we also show that there is an optimal 
localizing potential that minimizes the density oscilla- 
tions around the QHs. We also find that in systems with 
two QHs, in the Coulomb case, it is sufficient to have a 
single (5-function: This localizes a single QH in the GS, 
while the other is automatically at the antipodal point 
on the sphere, because of the Coulomb repulsion. 

We then study how states containing two and four QHs 
at 5/2 evolve as the interaction is varied between Vsf, and 
Vc, and show that there is AC between the two interac- 
tion limits, as in the previous work.l^l 

For two QHs, as the interaction is varied from Vsf, 
to Vc, we find that the lowest energy state does not 
mix with the higher lying ones. For systems containing 
four QHs, the two lowest-lying states forming a degener- 
ate doublet in the V^b limit'^^ ■^^ remain the two lowest- 
energy states even in the limit of pure Coulomb inter- 
action Vc- Thus, the Majorana fermion associated with 
the doublet'^^ appears to survive in the Coulomb limit. 

To better study the MR limit we then repeat the inves- 
tigations for four QHs using a different pinning potential: 
Also with this approach we obtain AC for the lowest-lying 
doublet between the two interaction limits, confirming 
the above results. Finally we present some preliminary 
results of braiding with QHs, showing that exchanging 
the position of two QHs, in the presence of two others 
fixed QHs, the system goes from one of the MR doublet 
states to the other, showing non-Abelian braiding statis- 
tics. Furthermore, we investigate the fusion of two QHs, 
getting an estimate for the MR doublet splitting.— 



II. THE SYSTEM 

In the following we consider fully spin-polarized elec- 
tron systems on the surface of a sphere:'*" By exact diag- 
onalization we obtain their low-lying states and energy 
spectra. In the spherical geometry for a half filled Lan- 
dau level (LL) the number of electrons TV and the number 
of fiux quanta are related by = 2N — S, where 
the shift S (Ref. 31) is a topological quantum number 



that depends on the particular FQH state. 

For our investigations of the FQH at filling factor 
1/ — 5/2 we study a half filled second LL, while the two 
lower- lying filled levels are considered as inert. For such 
a system the GS is at shift S = 3 (as obtained by exact 
diagonalization^i^, the same value as for the MR state.l^ 
We can insert Nqh e/4-charged QHs into the system by 
adding \Nqh supplementary flux quanta, i.e.. 



N^{N, Nqh) = 2N-'d + Nqh/2. 



(1) 



In the spirit of our previous worlP we consider inter- 
actions of the form 

V={l-x)V^P+ xV^"^ + Vp,n (x) , (2) 

where < a; < 1, interpolating between the Coulomb 
interaction Vc and the three-body interaction Vsb, whose 
exact GS is the MR wave function,!^ 
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^^b = ^ ^ S^,k {A,d{z - j)Ald{z - fc)} , (3) 

i<j<k 

when projected onto the lowest (spin-polarized) LL. Here 
Sijk expresses the full symmetrization over the permuta- 
tions within the triplet (ijk) and S{i — j) is the ^-function 
in the separation of particles i and j. The constant A is 
chosen as in the earlier work,!^ such that the gap for V^i, 
is approximately the same as that for Vc (note that V^b 
induces an extensive energy for all GS at any i' > 1/2). 

The superscripts (0) and (1) in Eq. Q indicate that 
the Coulomb energy V^^^ is evaluated in the n ~ 1 LL, 
as required for the experimentally realized v = b/2 state. 



while the three-body interaction Vgj^ 



(0) 



is evaluated in the 

n = LL, since the MR state is the GS of V^b in the 
lowest LL. 

The pinning term Vpin{x) localizes the QHs at fixed 
positions on the surface of the sphere. In our first 



calculations (in sections HI and IV A) we consider a 
parametrization of the pinning potential of the form 



VpUx) = {l-x)-[{l-x)V^^l + xV^ 



(4) 



where the pinning function Vpi„ is still to be chosen, 
such that it localizes exactly an e/4-charged QH. This 
ypin{x) couples to the electron density, interpolating like 
the electron interaction between LL n = 1 and n = 
with weights (1 — x) and x, respectively. The overall fac- 
tor [1 — x) makes the localization term vanish in the MR 
limit as x — > 1, just like the Coulomb term, both helping 
to separate the energies of states with charge e/4 QHs. 

However, in the MR limit {x — 1), the e/2-charged 
"double QH" has a zero in the electron density and thus 
zero pinning energy, while the e/4-charged QH shows a 
finite local minimum of the density and thus finite pin- 
ning energy. To separate e/4-charged QHs requires ei- 
ther an additional repulsive interaction between QHs or 
a special pinning potential coupling to the monopole and 
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Figure 1. (Color online) Polar angle dependence of electron density and missing charge for v — 5/2 systems with two QHs 
localized at the poles of the Haldane sphere, for Coulomb interaction. In (a) the QHs are localized by a (5-function potential Vpin 
[setting a; = in Eq. Q], for two different pinning strengths; in (b) by a smooth pinning potential, whose radial dependence 
is shown in (c). (d) Energy spectrum of a system with two QHs, without pinning, as function of the total angular momentum 
L, and comparison with the Coulomb energy for S function and smooth pinning [cf. the text for more details; results are for 
N = 20, = 38 in (a) and for N = 18, = 34 in (b),(c) and (d)]. 



quadrupole moment of the QHs (cf. problems discussed 
in Toke et alJ^. To avoid this difficulty, in this first se- 
ries of calculations we never reach exactly the MR point, 
but approach it closely enough to identify the lowest- 
energy states. In Section [IVB| we will then use a different 
localization potential'^'^ to circumvent this problem and 
study the full AC up to the MR point. 

All energies are measured in the usual units e'^/eio, 
where Iq = y/hc/eB is the magnetic length, and taking 
the energy of the lowest state as reference. 



III. LOCALIZATION OF QUASIHOLES 

We start by studying the QH localization in the pure 
Coulomb limit, taking for Vpin a (^-function pinning po- 
tential: Choosing a suitable pinning strength we see that 
it is actually possible to obtain separated e/4-charged 
QHs, in contrast to the results of Toke et al}^ Figure 
1(a) shows the polar angle dependence of the density (up- 
per curves) and of the integrated missing charge (lower 
curves) for two different pinning strengths, for a system 
containing two QHs with localizing 5 functions, one at 
the north and the other at the south pole of the sphere. 
We observe that the two QHs are localized at the poles. 
However they are large in extent and are accompanied 
by significant and slowly decaying density oscillations (cf. 
Nuebler et al^. The missing charge e/4 is reached (for 
the stronger pinning potential) at i? « 4£o and at the 
equator as required by symmetry. 

As the system is compressible at the position of a QH, 
the shape of the pinning potential will influence its struc- 



ture. We therefore study how the pinning potential can 
be tuned to minimize the size of the QH and to suppress 
the density oscillations around it. With the choice of 
shape, depicted in Fig. 1(c), we obtain an electron den- 
sity without oscillations, as illustrated in Fig. 1(b): The 
smooth potential allows to properly separate the two QHs 
and to reduce their radius to a minimum of ^ 'Mq. 

How is this possible? Figure 1(d) shows a comparison 
of the interaction energies, not including the energy due 
to the pinning forces, for (5-function pinning (upper line) 
and for smooth pinning (lower line). The circles depict 
the spectrum of a system with two free QHs, i.e., without 
pinning force, as function of the total angular momentum 
L. We see that the use of a smooth pinning potential 
allows to obtain a smaller interaction energy, admixing 
only the lowest-lying L eigenstates, while the 5- function 
potential admixes also higher- lying states, giving rise to 
the large density oscillations. We also note that a "better" 
localization, with a pinning potential that further lowers 
the interaction energy, is unlikely to be possible, because 
a minimal number of low-energy L eigenstates is needed 
to obtain the desired QH localization. 

We note that the energy of two-QHs systems is affected 
by how well the density oscillations fit with the separa- 
tion distance of the QHs, such that the energy will show 
oscillations as function of the QH separation. The en- 
ergy oscillations observed in the variational calculation 
of the Majorana fermion energjEH may to some extent be 
caused by such commensuration effects. An optimized 
localization potential might help to reduce them, hope- 
fully allowing a more accurate calculation of the coher- 
ence length associated with the Majorana fermion. 
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Figure 3. (Color online) Low-lying energy spectrum as func- 
tion of the interaction parameter x for a system of A'^ = 12 
electrons with two QHs localized at the poles by 5-function 
pinning potentials; the lowest-energy state is the reference. 



the state at — N/2 — 8 which corresponds to a "dou- 
ble QH" with charge e/2 localized at the south pole, but 
no charge deficiency at the north pole. At i' = 5/2, for 
Coulomb interaction, this "double QH" state has quite 
large interaction energy and, except perhaps in the limit 
of very strong pinning, is unlikely to become the GS. 



Figure 2. (Color online) (a) Energy spectrum of a two-QH 
system as a function of the angular momentum component Lz 
when S functions are inserted at both poles (for two different 
pinning strengths, — 16, A^^ — 30); the lowest-energy state 
is the reference, (b) The same for a single 5 function at the 
south pole. 



In the following, we use (5-function pinning potentials. 
We first study what happens if only a single (5-function 
is inserted in a system of iV = 16 electrons with two 
QHs, and whether they are still separated or if a sin- 
gle e/2-charged "double quasihole" is formed. Figure 
2(a) shows the energy spectrum as function of the an- 
gular momentum component when two S functions 
are inserted at the poles, for weak and strong pinning - 
cf. Fig. 1(a). In both cases the lowest-energy state has 
Lz ~ and describes two separated QHs localized at the 
poles. Figure 2(b) shows the same, but for a single S 
function at the south pole. The lowest energy state has 
also Lz = 0, indicating that the QHs are still at max- 
imal separation, which minimizes their Coulomb repul- 
sion: The QH pinned at the south pole repels the other 
to the north pole of the sphere. For increasing pinning 
strength some states with are reduced in energy. 

With strong enough pinning, a state with ^ 0, in 
which QHs would cease to be maximally separated, may 
become the lowest energy state. Of particular interest is 



IV. ADIABATIC CONTINUITY 
A. Quasiholes localization with <5-function pinning 

We now turn to the study of AC for systems containing 
localized QHs. We begin with two QHs, pinned at the 
poles by S functions, and we vary the particle interaction 
as in Eq. ([2| from a; = (the Coulomb case) to near 
X — 1 (the MR limit). Figure 3 shows the evolution 
of the excitation energies as function of the parameter 
x for TV = 12 electrons and = 22 flux units; other 
system sizes give similar results. The horizontal line at 
zero energy represents the GS energy: Clearly none of the 
lines for the excited states come close to it and there is no 
mixing between them. Thus there is no phase transition 
in going from the Coulomb case to the MR limit, exactly 
as observed for systems without QHs.*" 

Next we investigate a system of = 14 electrons 
and = 27 flux quanta that contains Nqh = 4 lo- 
calized QHs.l^ This is particularly interesting because in 
the MR limit the GS is a degenerate doublet,!^ which 
is associated with a Majorana fermion^S- and w ith the 
non-Abelian braiding statistics of the QHs|2S123 For any 
four-QH configuration there are two linearly independent 
wave functions that describe it; braiding QHs around 
each other induces a linear transformation in the degen- 
erate subspace of the doublet. Does this doublet survive 
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Figure 4. (Color online) Low-lying energy spectrum as function of the interaction parameter x for systems with four QHs 
localized with 5-function pinning potentials: (a) At the corners of a square on a great circle and (b) at the vertices of a 
tetrahedron; the curves labeled Pfaffian give the variational estimate of the doublet splitting obtained using the GS doublet 
for the MR limit, (c) Overlaps of the two doublet-states for Coulomb interaction with the corresponding states for interaction 
parameter x. (d) Low-lying energy spectrum for Coulomb interaction, as function of the pinning strength, for four QHs in 
tetrahedral position. (The lowest-energy state is the reference; = 14, A^^ = 27.) 



in the Coulomb limit? 

Figures 4(a) and 4(b) show results for four QHs local- 
ized at the corners of a square on a great circle and at 
the vertices of a tetrahedron, respectively. In both cases 
we see that the two lowest-energy states close to the MR 
limit, i.e., the GS at energy and the first excited state, 
remain the lowest-energy states over the entire parame- 
ter interval, even at the Coulomb point x — 0. During 
the whole evolution there is neither crossing nor mixing 
with higher-lying states, indicating the absence of a phase 
transition. Note that in the MR case, at a; = 1, all the 
considered states have zero energy, because the pinning 
potential V^m(x) vanishes at this point. 

In the tetrahedral case, right in the Coulomb limit, 
a higher-energy state comes down in energy and nearly 
reaches the upper state of the doublet. But Fig. 4(d) 
shows that this is actually not a problem: The Coulomb 
spectrum is plotted as function of the strength of the 
pinning potential and it is evident that the near degen- 
eracy is easily lifted by increasing the pinning strength. 
We also note that the two states of the doublet are not 
degenerate in the presence of a Coulomb interaction but 
we expect the splitting to vanish in the thermodynamic 
hmit.EH 

The curves denoted by "Pfafhan" in Figs. 4(a) and 4(b) 
show the variational results for the splitting of the dou- 
blet, computed using the states at x = 0.95 (compare 
with Baraban et al.-'^): They overestimate the splitting 
in the Coulomb limit, in the square geometry by a factor 
of 2.5 and in the tetrahedral geometry by ~ 70%. 

Finally, Fig. 4(c) shows the evolution of the overlaps 
(V'i(O)IV'i(a;)) of the two states ^i{x) {i — 0,1) of the 
doublet as the interaction varies from Coulomb Vc at 
a: = to the vicinity of the three-body limit V^b &t x — 
0.95. We see that there is no sign of an abrupt drop of 
the overlaps which could signal a phase transition, and 
the values reached near the MR limit are reasonably large 
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Figure 5. (Color online) Low-lying energy spectrum as func- 
tion of the interaction parameter x for systems with four QHs 
localized at the vertices of a tetrahedron using the pinning 
potential in Eq. ([5|: (a) = 12, 7V^ = 23, g = -e/4, (b) 
N — 14, A^^ = 27, g = — 0.19e (the lowest-energy state in the 
MR limit is the reference). 
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Figure 6. (Color online) Particle density on the sphere surface for the two states of the MR doublet: (a),(b) in the MR limit 
and (d),(e) in the Coulomb case; (c) and (f) show the density difference between the two states (A'^ — 14, A^^ = 27). 



for a system of this size (TV = 14, = 27). 

All these results provide evidence in favor of adiabatic 
continuity from the MR to the Coulomb limit for the two 
lowest-energy states in systems with four localized QHs, 
suggesting that the non-Abelian doublet, and the associ- 
ated Majorana fermion, can be expected in the limit of 
pure Coulomb interaction. 



B. Quasiholes localization with "STM tip" pinning 

In the AC investigations described so far, the pinning 
potential Vpin{x) [in Eq. (|4|] was not able to localize QHs 
in the MR limit (as discussed in section |ll|, but we ap- 
proached this limit near enough to recognize the lowest- 
lying states. One can pose the question if the permanent 
presence of the Coulomb interaction (that vanish as 1 — a; 
in the MR limit a; — > 1) could influence our results, se- 
lecting exactly the "right states" from the zero energy 
states at the MR point. To strenghten our evidences and 
to show a complete AC between Coulomb and MR limit, 
we repeat here the calculations for systems containing 
four QHs, using a different pinning potential Vpin{x) that 
allows to reach both limits: 



VpM = {I - x)Vj'^ + xV.^^'l 



(5) 



where Vq, interpolated between the n — I and the n ~ 
LL like the electron interaction, is the Coulomb potential 
of a pointlike object, e.g., a STM tip, with a (negative) 
charge q, positioned on the surface of the FQH sample. '^'^ 
This potential repels the electrons from its center and, 
if the charge q is chosen in an appropriate way, namely. 



in an interval near — e/4, it can only localize e/4-charged 
QHs, thus avoiding the problem of the localization of 
e/2-charged "double QHs" in the MR limit. This idea 
was first used by Prodan and Haldan^^ in their inves- 
tigation of the (non-Abelian) braiding properties of MR 
QHs, although projecting the pinning potential on the 
zero energy space of their Hamiltonian. 

Figure 5 shows the results for the low-lying energy 
spectrum, as the interaction is varied from the Coulomb 
{x = 0) to the MR limit (a; = 1), for systems with four 
QHs localized at the corners of a tetrahedron by four 
such pinning potentials. The results in Fig. 5(a) are for 

= 12 electrons and a localizing charge q = —e/4, while 
in Fig. 5(b) for = 14 electronpl and q = -0.19e: 
These charges were chosen in order to obtain that the 
two states of the MR doublet are exactly the lowest-lying 
states in the MR limit. 

In both cases we see that (once q is fixed) the two 
lowest energy states remain the lowest over the whole 
parameter interval (for A^ = 12 with a level crossing). 
The third energy level, corresponding to an excitation in 
the pinning potential, is exactly threefold degenerate (for 
symmetry reasons) and can thus be easily distinguished 
from the two lowest states forming the MR doublet. 

We note that this doublet has a finite energy splitting, 
not only for pure Coulomb interaction, but also in the 
MR limit: This is caused by the pinning potential ^ 
that mixes in states with nonzero energy. We also find 
that the energy separation between the upper state of the 
doublet and the third (threefold degenerate) state can be 
small, particularly in the MR limit, and depends on the 
charge q used in the pinning potential. For example for 
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N = 14 electrons, using a localizing charge q = — e/4. 
the pinning potential is too strong and at the MR point 
the threefold degenerate state is slightly lower in energy 
than the upper state of the doublet; to obtain the right 
order we chose a weaker potential, with q — — 0.19e. 

The AC between the Coulomb and the MR limit for the 
lowest energy doublet shown in Fig. 5 thus supports the 
evidences of the previous section in favor of a "survival" 
of the Majorana fermion till the Coulomb point. 

In Fig. 6 we also show the particle density on the sur- 
face of the sphere for the two states of the MR doublet in 
the two limiting cases: (a), (b) for the MR and (d), (e) 
for the Coulomb limit. The two states of the doublet are 
indeed very similar, as also shown in (c) and (f), where 
the density differences are plotted. Note that the density 
differences in the MR case are approximately twice as 
large as the differences in the Coulomb limit. 



V. QUASIHOLES BRAIDING AND FUSION 

In this last section we show some preliminary results 
of QH braiding and fusion in systems containing four lo- 
calized QHs. We perform this investigation for electron 
interactions near to the Coulomb limit jbut slightly in- 
creasing vi, the Haldane pseudopotentiaP^lthat describes 
the interaction between two particles having a relative 
angular momentum Ih. In our previous work^ we showed 
that, under such interaction modifications, the system re- 
mains in the MR phase, even improving the overlap be- 
tween the exact GS and the Pfafhan wave function and 
simultaneously with an enhancement of the gap. For the 
QH localization we use the pinning potential introduced 
in the previous section [Eq. ([5|], setting x — 0, that is 
evaluating it fully in the n = 1 LL, and chosing a suitable 
value for the localizing charge q, as explained below. 

Investigating QH braiding in systems with four QHs, 
the non-Abelian nature of the MR doublet states should 
become manifest: Interchanging the positions of two QHs 
by stepwise changing the location of their pinning po- 
tentials (and keeping the other two QHs fixed), the two 
states of the MR doublet should transform into each 
other, in particular, switching their position in the en- 
ergy spectrum. For this to happen, during the QHs in- 
terchange process, an odd number of level crossings be- 
tween them is needed. In the following calculations, we 
engineer an exact energy degeneracy for the two states of 
the doublet at the midpoint of the QHs interchange pro- 
cess, by choosing suited values for the localizing charge q 
and the first Haldane pseudotential vi. Then, performing 
the interchange, we investigate if this is indeed a crossing 
point and whether it is the only one or if others arise. 

We chose as starting (and ending) configuration four 
QHs at the corners of a rectangle: Two of them on the 
upper half-sphere, at the same polar angle Oi^2 = = 
02 = 54.736° and at the opposite azimuthal angles ipi ~ 
90° , (p2 — 270° ; the other two on the lower half-sphere at 
the polar angle 63,4 = 63 = O4 = 125.264° and azimuthal 
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Figure 7. (Color online) Braiding of two QHs (in the pres- 
ence of two other fixed QHs): low- lying energy spectrum as 
function of the braiding angle (p, defined through Eq. ([6|. 
The electron interaction is near to the Coulomb limit (with a 
slightly increased Haldane pseudopotential Vi) and the QHs 
are localized by the pinning potential of Eq. ([5]), setting x — 0. 
(a) TV = 12, q = -0.1715e, vi = 1.05uf°^, (b) iV = 14, 
q = -0.19e, vi = 1.075wf°"', where is the Coulomb 

value of the first Haldane pseudopotential. 

angles ip3 = 90° and (^4 = 270°. We keep QHs 3 and 
4 fixed at their locations and, by stepwise changing the 
positions of the pinning potentials, we rotate both QH 1 
and 2 around the vertical axis through the poles, 

ipiiip) = 90° + if, ip2i(p) = (270° + if) mod 360°, (6) 

keeping the polar angle 9i 2 — 54.736° constant, until 
they exchange their original positions: 0° < (p < 180°. 
This is a 180°-periodic process; we call ip the "braiding 
angle". 

The polar angles 612 and 63 4 are chosen in such a 
way that at the midpoint of the rotation process, that 
is at braiding angle ip = 90°, the four QHs are in a 
tetrahedral configuration. The localizing charge q and 
Haldane pseudopotential vi (for the whole braiding pro- 
cess) were previously chosen, such that at this point 
the two states of the MR doublet are degenerate in en- 
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ergy: q = -0.1715e, vi = 1.05i;f°"' for TV = 12 and 
q = -0.19e, vi = l.OTS'yf for = 14 electrons, where 
yCoui jg |.]^g Coulomb value of the pseudopotential vi. 
These are "reasonable" values for q and vi: The local- 
ized QHs have charge e/4 and the electron interaction is 
modified such that the system is still in the MR phase. 

Figure 7 shows the evolution of the low lying energy 
spectrum as function of the braiding angle ip for (a) 
= 12 and (b) = 14 electrons. In both cases we 
observe that the two states of the doublet remain the 
lowest-energy states during the whole braiding process, 
never mixing with higher-lying states or with each other 
(except at the degeneracy point). Indeed the overlaps 
between the corresponding states of the doublet for suc- 
cessive braiding steps are consistently high: In the range 
0.94-0.97 for TV = 12 (5° steps) and 0.97-0.99 for = 14 
(3° steps). 

From the plots it is evident that the degeneracy point 
a,t (fi — 90° is indeed also a crossing point for the two 
doublet states: If one follows the lines with a continuous 
derivative, at the degeneracy point one goes from the 
lower to the upper state of the doublet (and vice versa). 
No other crossing or mixing point arises and thus the QH 
braiding causes an interchange in the MR doublet states: 
If one starts aX (p — 0° in the lower-energy state, one ends 
a.t (p — 180° in the upper-energy state, as expected for 
non-Abelian braiding statistics. Two such braidings are 
needed to come back to the initial state. 

We note that we can get an exact degeneracy of the 
MR doublet only when the four QHs are in the tetra- 
hedral configuration. Slightly modifying the geometry 
of the configuration, a small gap opens; however, if the 
speed of the braiding is sufficiently high, the states of the 
MR doublet would still interchange at this point. This 
situation is analogous to that studied by Thouless and 
Gefen,"^^ concerning the crossing between the lowest-lying 
states as function of the magnetic fiux, for the quantum 
Hall effect at fractional fillings, showing that these cross- 
ings are essential in order to get a fractional charge. On 
the other hand if the QHs are moved in a strictly adia- 
batic way, the system always remains in the lowest-energy 
state, without crossing. We also wish to emphasize that 
the doublet degeneracy is very different from the three- 
fold degeneracy of the first excited state: The former 
is obtained only by fine tuning of the parameters q and 
vi, the latter is a purely geometric degeneracy, resulting 
from the symmetric QH configuration and independent 
from interaction and localization parameters. 

The fact that the MR doublet states are not degener- 
ate (except by tuning at ip = 90°) in our small system 
diagonalizations helps us to follow their evolution during 
the braiding (through the evolution of the energy lev- 
els) and thus to recognize the crossing point which leads 
to the interchange of states. However we expect that in 
the thermodynamic limit the two states become exactly 
degenerate during the whole braiding process, when the 
QHs are sufficiently far apart from each other: On one 
side this would make the fine tuning of the parameters q 
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Figure 8. (Color online) Fusion of two QHs at the north pole 
(in the presence of two other fixed QHs): Low-lying energy 
spectrum for A'' = 14 electrons as function of Si, 2 = d\ = 82, 
the polar angle of both the moving QH 1 and 2. The elec- 
tron interaction is near to the Coulomb limit (with a slightly 
increased Haldane pseudopotential v\ = l.OTSwp""') and the 
QHs are localized by the pinning potential of Eq. (|5|, setting 
a; = and q = — 0.19e (t;f is the Coulomb value of the first 
Haldane pseudopotential) . 



and vi unnecessary, but on the other side it would make 
it impossible, using our present method, to follow the 
evolution of the states under braiding (however, it is still 
possible with a study of the monodromy matrix, as was 
done by Prodan and Haldan^^filjj nevertheless, the inter- 
change of the MR doublet states should survive also in 
this limit. 

Unfortunately we can get results for QH braiding only 
for small system sizes '^^ and from them we cannot extract 
useful information about the dependence of the energy 
splitting when approaching the thermodynamic limit. 
The average splitting of the MR doublet for A^ = 14 
is slightly smaller than for A'^ = 12, but this is mainly 
caused by the different shape of the energy oscillations 
as function of the braiding angle Lp in the two cases: For 
A^ = 14 the oscillations of the MR doublet states are al- 
most "in phase", while for A^ = 12 they are "out of phase", 
thus giving a larger average splitting. 

Finally we investigate the fusion of two QHs by bring- 
ing the pinning potentials that localize them close to each 
other, in the presence of two others fixed QHs. We start 
from the tetrahedral configuration described above and 
we let fuse the two QHs in the upper half-sphere, by tak- 
ing them to the north pole, that is, by shrinking the polar 
angle 6*12 ^ Oi = O2 from 54.736° to 5° at fixed azimuthal 
angles. Figure 8 shows the low- lying energy spectrum as 
function of the polar angle 6*1 2 during this process, for 
A^ = 14 electrons. We see that, after that the degeneracy 
is lifted, the evolution of the two lowest states is similar, 
with a relatively constant splitting between them. This 
splitting (or at least a part of it) could come from the in- 
trinsic splitting of the MR doublet as two QHs are taken 
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together (its order of magnitude is the same as that ob- 
tained by Baraban et al.^^). However, at approximately 
9i 2 = 37°, the upper state of the MR doublet comes very 
near to the next higher-lying state and strongly mixes 
with it'^^ (but without crossing) , possibly losing a part of 
its character. Thus it is not fully clear what happens at 
the MR doublet, and further investigations are needed. 

VI. CONCLUSION 

Our results provide evidence in favor of adiabatic conti- 
nuity from the Moore-Read to the Coulomb limit for the 
two lowest-energy states in systems with four localized 



QHs and thus the non-Abelian doublet and associated 
Majorana fermion can be expected in the limit of pure 
Coulomb interaction. Forces breaking particle-hole sym- 
metry like V^t are actually present due to LL mixin^^^J^ 
and will favor either the Pfafhait^or the Antipfafhan.'^l^ 
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